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ABSTRACT 

It has been recently realized (Rephaeli 1995) that the relativistic corrections 
to the spectral distortions of the cosmic microwave background (CMB) measured 
in the direction of clusters of galaxies containing hot gas are significant and 
should be detectable with the forthcoming experiments. In the present paper 
we calculate the correction terms that are proportional to V r /c x kT e /m e c 2 and 
(V/c) 2 to the standard formulae describing the spectral distortions caused by 
the bulk motion of the free electrons (kinematic effect) and due to the presence 
of the hot gas (thermal effect) for the case of a cluster having a peculiar velocity 
V (V T is its radial component). The results of our analytical calculations are 
confirmed by Monte-Carlo simulations (Sazonov & Sunyaev 1998). 

Subject headings: Cosmology: theory — cosmic microwave background radiation 

- galaxies: clusters: thermal and kinematic Sunyaev-Zel'dovich effects - 
plasmas: Compton scattering 
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1. Introduction 

Thomson scattering of the cosmic microwave background (CMB) radiation by hot 
electrons in the intergalactic gas in clusters of galaxies modifies the spectrum of the CMB 
(Sunyaev and Zel'dovich 1972). Zel'dovich and Sunyaev (1969), basing on the Kompaneets 
equation (1957), derived a simple formula describing the spectral form of the distortion, 
which is proportional to the parameter y = (kT e /m e c 2 )r, where r is the Thomson optical 
depth along the line of sight. The effect has now been observed from a number of clusters 
of galaxies (see Birkinshaw 1998 for review). 

Recently, interest to this effect has been reactivated in view of the perspectives 
of accurate measurement of the CMB distortions in a number of experiments, both 
ground-based and on balloons, by the MAP spacecraft and especially by the Planck 
Surveyor mission scheduled to be flown in the middle of the next decade. These activities 
were motivated by the fact that the gas temperature is so high in the clusters of galaxies 
(ranging between 3 and 17 keV, Tucker et al. 1998) that the scattering electrons have 
thermal velocities of the order of 0.1 - 0.3 c, so one has to include into consideration the 
relativistic corrections to obtain an accurate result. Rephaeli (1995), basing on extensive 
previous work (Wright 1979, Fabbri 1981, Taylor & Wright 1989, Loeb et al. 1991) has 
demonstrated by means of numerical calculations the relevance of the relativistic corrections 
for the future experiments. Stebbins (1997), Itoh et al. (1998), and Challinor & Lasenby 
(1998) used a Fokker-Planck approximation of the relativistic photon kinetic equation to 
obtain corrections, written as series in powers of kT e /m e c 2 , to the standard nonrelativistic 
solution. These results have proved to be in excellent agreement with those of Rephaeli, 
demonstrating the applicability of the diffusion approximation to the problem at hand, 
despite the small optical depths of the clusters of galaxies (r ~ 0.01). 

A gas cloud moving rapidly relative to the CMB along the observer's line of sight 
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must significantly modify the spectrum of the CMB in addition to the thermal effect. 
The change in the brightness temperature caused by this "kinematic" effect is to first 
order simply proportional to the radial component of the cluster velocity ~ (y r /c)r 
(Sunyaev & Zel'dovich 1980). The effect should be detectable in the future, and will enable 
measurement of cluster peculiar velocities, with significant implications for studies of the 
large-scale structure of the universe. It is obvious that corrections similar to those found 
for the thermal effect must exist and should be taken into account for the kinematic effect, 
if one wants to find the correct solution for the case of a moving cluster. In this paper 
we calculate the next-order changes in the spectrum of the CMB related to the cluster 
peculiar velocity by solving the photon kinetic equation. We have obtained simple formulae 
giving the correction terms of the orders of V T /c x kT e /m e c 2 and (V r /c) 2 . Our method is 
similar to that used by Psaltis & Lamb (1997) who considered the more general problem of 
comptonization in a moving media. The solution these authors have obtained, although 
applicable to many astrophysical situations, does not contain the 0(V r /c x kT e /m e c 2 ) term, 
because this term is third-order in electron velocity, whereas their solution is accurate only 
to second order in it. We aslo confirm the existence of the term of order (kT e /m e c 2 ) 2 found 
earlier using techniques different from ours (Rephaeli 1995, Stebbins 1997, Itoh et al. 1998, 
and Challinor & Lasenby 1998). Earlier we have found all the correction terms mentioned 
above numerically using Monte-Carlo simulations (Sazonov & Sunyaev 1998). 



2. Scattering of the CMB by a directed beam of electrons 

From the point of view of the observer each electron in the intergalactic gas scatters 
the CMB photons independently. We can therefore first consider the problem of scattering 
of the CMB by a directed beam of monoenergetic electrons having a density N e and moving 
at velocity v. Once we have found a solution to this simplified problem, we will be able to 
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consider the more general problem of scattering of the CMB by a cloud of thermal electrons 
having a peculiar motion, which corresponds to the real situation of a cluster of galaxies, 
by simply averaging the result obtained for the directed beam over a drifting Maxwellian 
distribution of electron velocities. We will use two coordinate frames. Quantities with 
subscript refer to the system where the electrons are at rest, while quantities without 
subscript refer to the frame that is fixed to the CMB, hereafter referred to as the laboratory 
frame. It is easy to show that our results are valid for a cluster at any redshift (see Sazonov 
& Sunyaev 1998). In the laboratory frame the initial occupation number in the photon 
phase space is planckian with a temperature Tcmb: n = l/(e x — 1), where x = hu //cTcmb- 
The corresponding spectral intensity is 

_ 2(fcT CMB ) 3 x 3 

{hcf e*-l [) 

The occupation number is invariant with respect to the Lorentz transformations of the 
frequency and direction of motion of a photon 



u ° = n ^ n = 7(1 - PfJ>) and /i = ^ , (2) 

where (5 — v/c, 7 = (1 — /3 2 ) -1 / 2 , li = cos 9, and 9 is the angle between v and the photon 
velocity (Landau & Lifshitz 1975). Therefore, in the electron rest frame the occupation 
number depends on the photon incident direction 



Ignoring induced scattering and the change in the photon frequency in each scattering 
due to the recoil effect, we can write down the photon kinetic equation in the electron rest 
frame as follows: 
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-^-n (/i , Vq) = cN c0 J J^MaC v ) - n (/i , fo)] dO , (4) 

Integrating the Thomson differential cross-section over the azimuthal angle (Chandrasekhar 
1950) one derives 



-£-no(»o, ^o) = 3c ^ c ° aT f (3 + 3/ig/iQ 2 - nl - ^) M/4 u ) - noOo, ^o)] d/j,' (5) 

Clio J-DTT "'-I 

We now evaluate the collision integral in equation (|5[) by expanding no up to forth 
order in (3. The integral is then easily taken, yielding: 
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Our next step is to calculate the corresponding scattering rate as measured in the 
laboratory frame. This can be done by making use of the Lorentz-invariance property of 
the photon occupation number. It is easily shown (see e.g. Peebles 1971) that 



c/n(/i, v) 



dn (fj, , u ) 



(7) 



dt 7(l + /3yU ) dt 

Using relations (0), (|7]) and substituting N e0 = N e /^ (due to the Lorentz-transformation 
of volume) for the electron density, we derive from equation (Qj 
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dn(fi, v) 



cN e a T — 



xe 



n(/j,, v) dt c e x ' — 1 

9 + 13/i 2 



-1 -/i 2 + 



-11/i 2 
20 



" 31 + 11/i 2 
20 



+ 



120 



-(2F 2 + G 2 



+ P 4 



■2 " ■>■>!• 9 + 66/x 2 + 13// 2 ^ ^ 2 



20 120 
3 + 33/x 2 + 28/x 4 
+ 420 



F(F 2 + 2G 2 



where F = xcoth (x/2), and G = x/ sinh(x/2). 



Compton scattering must save the total number of photons. We have verified that all 
the (3 terms in equation (|8|) indeed vanish after the integration over photon direction and 
frequency: d/dt J nv 2 dv d[i = 0. 

Another known integral property of the process of compton scattering is the energy 
exchange rate between an electron and an isotropical radiation field. The radiation energy 
density e r (see e.g. Pozdnyakov et al. 1983) should increase with time as 



de x 
~~dl 



-cN e a T e r (-f 2 - 1) 



(9) 



Taking an integral d/dt J nv 3 dud[i we find in our case 



§ = ;W T£r (/? 2 + n 

dt 3 

which is identical to dependence @ to forth order in (3. 



(10) 



3. Scattering of the CMB by the hot gas in a moving cluster 

Consider now a cluster of galaxies moving with a peculiar velocity V at an angle 9 
(fi = cos 6) relative to the vector drawn from the cluster to the observer. The cluster 
contains hot gas, and the distribution of the electrons in the cluster rest frame is assumed 
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to be relativistic Maxwellian with a temperature T e : dN e = Aexp [— E (p Q )m e c 2 / kT 2 } dpo, 
where po is the electron momentum, E is the electron energy, and A is the normalization 
constant. The corresponding distribution in the laboratory frame is obtained via the 
Lorentz-transformation of p 



Px = liPxO + ^ E o)', Py = PyO] Pz = PzO, (H) 

where 7 = (1 — V 2 / c 2 )^ 1 / 2 , and axis X is drawn along the direction of the cluster peculiar 
motion (Landau & Lifshitz 1975). We can average equation ([8]) over the resulting electron 
velocity distribution expanded in powers of Vj c and T e to get 
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(12) 



Here we have replaced an integral over time / N e {r)a^cdt by an integral along the line of 
sight t = J N e (r)aTdr, where r < 1. A Monte-Carlo computation that takes into account 
only single-scattering events proves the validity of this transition (see a detailed discussion 
in Sazonov & Sunyaev 1998). 



The subsequent terms in equation ([12] ) physically correspond to increasing orders in j3 
in equation @, which is in the current problem a sum of the electron thermal and peculiar 
velocities. The first term (of order V/c) describes the kinematic effect. The second term, 
which is proportional to the second power of the thermal velocity, describes the thermal 
effect. The 0[(V/c) 2 ] term is the relativistic correction to the kinematic effect for a cloud 
of cold electrons. The "interference" term that is proportional to V/c x kT e /m e c 2 draws 
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from the term of order 



j3 3 in equation (^|). We have ascertained its existence by means of 



Monte-Carlo simulations (Sazonov & Sunyaev 1998). This term constitutes the leading 
relativistic correction to the kinematic effect and is the main subject of the present paper. 
Finally, the term of order (kT e /m e c 2 ) 2 (forth-order in (3) is the relativistic correction to the 
thermal effect found earlier using a Fokker-Planck approximation (Stebbins 1997, Itoh et 
al. 1998, Challinor and Lasenby 1998). 



The distortion of the spectral intensity of the CMB is related to the corresponding 
change in the photon occupation number by the equation 



In Fig. 1 we plot this distortion at kT e /m e c 2 = 0.02, V/c = 0.01 (the large value for 
the peculiar velocity has been chosen for illustration purposes), for two opposite directions 
of the cluster motion \i = 1 and \x = — 1. One can see that the contribution from the newly 
found term 0(V/c x kT e /m e c 2 ) to the total effect is significant. This contribution reaches 
its maximum at x — 3.34, i.e. near the frequency x c ~ 3.83 at which the thermal effect 
vanishes and where measurements of the kinematic effect would seem most promising. 

In our previous paper (Sazonov & Sunyaev 1998) we have performed Monte-Carlo 
simulations to determine the spectral changes in the CMB for various sets of parameters. 
As evident from Fig. 1, which uses the result of that paper, the correctness of the analytical 
formula (|12|) is confirmed by the numerical calculations. 

Using equation (|T2"D we have obtained a simple approximation formula describing the 
position of the crossover frequency (Xq), i.e. the frequency at which the distortion of the 



4. Properties of the CMB Spectral Distortion 
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incident microwave spectrum is zero 



V r m e c 2 _ _ kT P „ „ V r 



X = 3.830 1 - 0.31 — -^— + 1.1 \ - 0.6— (14) 

\ c kT e m e c 2 c j 



The correction term of order kT e /m e c 2 in this formula was known before our work (Rephaeli 
1995, Itoh et al. 1998, Challinor & Lasenby 1998). The 0(V/c x m e c 2 /kT e ) term is due to 
the kinematic effect. 

We have obtained similar approximation formulae for the positions [X min and X max ) 
and values (J mm and J max) of the minimum and maximum of the spectral dependences 
shown in Fig. 1. They are as follows 



/ V r m e c 2 kT e V T \ 
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X max = 6.511 (l - 0.09^^ + 2.5^ - 0.5^1 , (17) 

\ C KJ- g Tfl^C C J 



J ™ = 3 - 390 ~wr^ T i 1 + a43 7 -m ~ 6 - 2 ^ + °- 6 7 J ' (18) 

Finally, one can calculate the excess in the CMB energy flux in the direction of the 
cluster (51 = J 5I V dv) 



SI = tI 
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(19) 
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where I = 6Tq MB /(47tc), and b is the Stefan-Boltzmann constant. 

One can see again that the relativistic correction of order V/cx kT e /m e c 2 is important. 
Equation fll9|) allows one to calculate the rate of energy exchange between the hot gas and 
the CMB. We immediately see that the 0(kT e /m e c 2 ) and 0[(kT e /m e c 2 ) 2 } terms are the two 
leading terms in the relativistic formula that gives the energy transfer rate averaged over a 
Maxwellian distribution of electron velocities: de r /dt = 4/3 cN e aT^r(l 2 — 1)- Integration of 
the term of order iV/c) 2 in formula ([19]) over the observing angle (/i) leads again to relation 
(|i~0l) for the energy transfer rate due to the bulk motion of the electrons, as one should have 
expected. 

We have also obtained approximation formulae giving separately the flux from the 
"negative" source, i.e. integrated over the frequency range (0, X ) and that from the 
positive source (integrated from X to oo), taking into account dependence (|H]) 
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Fig. 1. — (a) An example of the CMB spectral distortion (in units of 2(kT e f / (he) 2 ) due to 
the Sunyaev-Zel'dovich effect for the following set of cluster parameters: kT e /m e c 2 = 0.02 
(T e = 10.2 keV), V/c = 0.01 (V = 3000 km/s), and p = 1 (the cluster moves toward 
us). The solid line shows the cumulative effect, which was calculated by summing the 
analytical expressions for terms of different orders in kT e /m e c 2 and V/c as given in Itoh 
et al. (1998) and the present paper. The contributions from the following components are 
shown: 0(kT e /m e c 2 ) (dotted curve), 0[{kT e /m e c 2 )] 2 (short-dashed curve), 0(V/c) (long- 
dashed curve), 0[(V/c) 2 ] (dash-dotted curve), and 0(V/c x kT e /m e c 2 ) (long-short dashed 
curve). For comparison, the result of Monte-Carlo simulations is shown as a histogram, (b) 
Same as (a), but the cluster moves outwards from the observer (// = —1). For comparison, 
the result shown in Fig. 1 is repeated. 
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